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Abstract In this chapter we introduce a microscopic modelling of the surplus elec- 
trons on the plasma wall which complements the classical description of the plasma 
sheath. First we introduce a model for the electron surface layer to study the qua- 
sistationary electron distribution and the potential at an unbiased plasma wall. Then 
we calculate sticking coefficients and desorption times for electron trapping in the 
image states. Finally we study how surplus electrons affect light scattering and how 
charge signatures offer the possibility of a novel charge measurement for dust grains. 



1 Introduction 

When a macroscopic object is brought into contact with an ionised gas its surface is 
exposed to a fast influx of electrons. Due to the collection of electrons it acquires a 
negative charge. The wall electrons give rise to a repulsive Coulomb potential which 
reduces the further influx of electrons until — in quasi-stationarity — electron and ion 
flux onto the surface balance each other This coincides with the formation of the 
plasma sheath, an electron depleted region adjacent to the wall. 

The traditional modelling of the plasma wall focuses mainly on the plasma sheath 
which is a macroscopic manifestation of the charge transfer at the plasma boundary. 
Electrons and ions are usually assumed to annihilate instantly at the wall. Thus the 
wall is considered to be a perfect absorber for the charged species of the plasma. 
This boundary condition implicitly assumes that the charge transfer at the plasma 
wall occurs on time and length scales too short to be of significance for the devel- 
opment of the discharge. It is thus not necessary to track them 1 1]. Improvements to 
the perfect absorber condition characterise the wall by an electron-ion recombina- 
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Fig. 1 Three novel bounded plasmas where surface charges matter. Left: The particle charge con- 
trols interactions between dust grains forming a spherical dust ball. Middle: Surface charges (neg- 
ative charges in blue over neutral background in red) pin the location of filaments for a filamentary 
dielectric barrier discharge. Right: The plasma serves as a collector in a plasma transistor (repro- 
duced from (2l)- 



tion coefficient and secondary electron emission coefficients for various impacting 
species. Still, replacing the electron kinetics by an ad-hoc boundary condition does 
not allow for a description of the charge transfer at the surface. 

While this approach may be sufficient for large discharges with small surface to 
volume ratios it becomes questionable for the modelling of small discharges with 
large surface to volume ratios. Indeed, the nature and build-up of surface charges 
have received increasing attention in various novel set-ups of bounded plasmas (il- 
lustrated in Fig.[T]i. In dusty plasmas LI. 4j^ |, for instance, the charge soaked up by 
the dust particles controls the overall charge balance of the discharge IS] . Moreover, 
the particle charge controls the interaction between dust particles and between dust 
particles and external electromagnetic fields. The particle charge is thus a central 
quantity of interest which should be known as precisely as possible. In dielectric 
barrier discharges 17] [8] |9] \lU\ |TT] IT2 l or microplasmas |fT3] [141 surface charges 
determine the spatio-temporal structure of the discharge ifTSl . Surface charges, for 
instance, may pin the filaments in the filamentary mode of the dielectric barrier dis- 
charge and give the discharge a memory across several cycles. Lastly, in solid-state 
based microdischarges 1.16. 17,1. the miniaturisation of the discharge increases the 
surface to volume ratio to such an extent that the (biased) plasma wall becomes an 
integral part of the discharge. The charge transfer from the plasma into the solid and 
vice versa may even be used explicitly, as for instance in the plasma transistor |2|. 

While the most prominent effect of charge transfer at the plasma wall is the 
plasma sheath the microscopic processes responsible for it cannot be captured in 
any sheath model which is solely based on the long range Coulomb potential and 
classical mechanics. Inspired by Emeleus and Coulter ifTSl as well as Behnke and 
coworkers fl^ 'T^ who envisaged a surface plasma of ions and electrons coupled to 
the bulk plasma by phenomenological rate equations characterised by sticking coef- 
ficients, residence times, and recombination coefficients we consider the build-up, 
distribution and release of surface electrons from a surface physics perspective. A 
sketch of the wall bound-processes which cannot be described in a sheath model is 
given by Fig. |2] The description of the plasma boundary has thus to be augmented 
by an interface region — the electron surface layer (ESL) — which encompasses not 
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Fig. 2 Illustration of the plasma and its boundary. The bulk plasma (left) is characterised by quasi- 
neutrality. In the sheath ions outweigh electrons, quasi-neutrality does not hold and a potential 
gradient is formed. Note that the bulk as well as the sheath have macroscopic extent and can be 
described classically. The microscopic processes at the surface, for instance electron reflection, 
sticking or desorption as well as electron ion recombination occur on a length scale smaller than 
the shortest collision length of the plasma and have to be described by quantum mechanics. 

only the Coulomb potential but also short ranged surface potentials. This gives a 
framework for identifying electron trapping sites, determining the electron distri- 
bution and the potential across the interface in quasi-stationarity and calculating 
electron sticking coefficients. Moreover it serves us to discuss the optical properties 
of the electron adsorbate on a dust grain. 

In the following section we will introduce our model for the ESL for an unbi- 
ased plasma wall. In section|3]we will calculate sticking coefficients and desorption 
times, that are average residence times at the surface, for the case where electron 
trapping takes place in the image states in front of the surface. Finally in section|4] 
we apply our microscopic model to the electrons on a dust particle in order to study 
charge signatures in light scattering by a spherical dust particle and propose thereby 
a novel optical charge measurement. 



2 Electron surface layer 

After the initial charge-up is completed the plasma wall carries a negative charge. 
In this section we extend the modelling of the plasma boundary to the region inside 
the solid and calculate the plasma-induced modifications of the potential and charge 
distribution of the surface within a quasi-static model. Although knowing the poten- 
tial and charge distribution across the interface may not be of particular importance 
for present day technological plasmas, it is of fundamental interest from an interface 
physics point of view. In addition, considering the plasma wall as an integral part of 
the plasma sheath may become critical when the miruaturisation of solid-state-based 
plasma devices continues. 

The current modelling of the plasma wall focuses on the plasma sheath and gives 
the potential as well as the electron and ion densities across that region. Our model 
for the ESL extends this into the solid ll20l . Specifically, we calculate the potential as 
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Fig. 3 Qualitative sketch of the charge double layer fonned at a plasma wall. The electrons de- 
pleted from the plasma sheath are accumulated in the ESL. Note that the boundary between sheath 
and ESL is located in front of the crystallographic interface. 



well as the distribution of wall bound electrons. Moreover, we give the spatial posi- 
tion, width and chemical potential of the surplus electrons. For this we use a graded 
interface potential ll2Ti to interpolate between the sheath and the wall potential and 
distribute the surplus electrons making up the wall charge in this potential under the 
assumption that at quasi-stationarity they are thermalised with the wall t22ll . 

As depicted in Fig. |3] we consider an ideal, planar interface at z = with the di- 
electric occupying the half-space z < and the discharge occupying the half-space 
z > 0. At the moment we focus on the physical principles controlling the electronic 
properties of the plasma-wall interface. Chemical contamination and structural dam- 
age due to the burning gas discharge are discarded. In the model we propose the 
wall charge to be treated as an ESL, which is an interface specific electron distribu- 
tion thermalised with the solid stretching from the plasma sheath over the crystallo- 
graphic interface into the bulk of the dielectric. 

The boundary between the ESL and the plasma sheath at zo is located in front of 
the surface. It is the position where the attractive force due to the surface potential 
(^suif ™d the repulsive force due to the sheath potential ^sheath balance each other. 
Thus, zo is given by 

0surf(2o) + 0sheath(2o) =0. (1) 

For z < zo an electron is attracted to the surface and thus contained in the ESL while 
for z > Zo it is repelled back into the plasma. The position zo is an effective wall for 
plasma electrons and ions at which the flux balance between electron flux je and ion 
flux ji has to be fulfilled. In the following we use for simplicity the perfect absorber 
model je — ji. The effective wall position also gives the distance to the surface at 
which the description of the plasma sheath based on the long-range potential breaks 
down as the short-range surface potential becomes dominant. On the soUd side, for 
z < 0, the ESL is bounded because of the shallow potential well formed by the 
restoring force from the positive charge in the plasma sheath. 

Within this quasi-static model the electrons missing in the sheath are accumu- 
lated on the surface. For the consttTiction of our one-dimensional interface model 
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we need the total number per unit area of missing sheath electrons (that is, the total 
surface density of missing sheath electrons) as a function of the wall potential be- 
cause it is this number of electrons which can be distributed across the ESL. Hence, 
we require a model for the plasma sheath. 

For simpUcity, we use a collision-less sheath model [l], more realistic sheath 
models IT] |23] |M1 make no difference in principle. In the collision-less sheath 
electrons are thermalised, that is, the electron density tie ~ noexp{e(j) /ksTe), with 

the potential, «o ths plasma density and Tg the electron temperature. The ions 
enter the sheath with a directed velocity v,o and satisfy a source-free continu- 
ity equation, d(njV,)/dz = 0, implying n,v/ = nov';o, and an equation of motion 
Mv,dv,/dz = — ed0/dz, with n,- the ion density, and M the ion mass. The poten- 
tial satisfies Poisson's equation d^0/dz^ = —4ne{nj — rie). Thus, the governing 
equations for the collision-less plasma sheath are [Tl 

dv,- e d0 



and 

d^ 
dz 



2 = —AneiiQ 



exp — - 



(3) 



Introducing dimensionless variables Tj = —e(j)/{kBTe), £, — z/Xd and u = Vj/cg, 
where Xd ~ \/ kgTe / Aitnoe^ and c., = ^Jk^TelM, equations (|2]iand(|3} become 

uu' — rj' and Tj" = — — exp(— Tj). (4) 
u 

In the ESL model the plasma occupies not the whole half space z > but only 

the portion z> zq (see Fig.jsj. Integrating the first equation gives u = —^2r\ +Mq, 

where mq — v,o/Ci is the reduced velocity of ions entering the sheath, so that the 
second equation becomes 

r]" = - J^ -exp(-r?). (5) 
'277 + Mq 

Using the boundary condition that the potential and the field vanish far inside the 
plasma, that is, Tj and Tj' -H- for oo, Eq. (jsjl can be integrated once and we 
obtain 



I?' = -\j -2uQ^J2rl+ul + 2ey.p[-rl)+2uQ^Jul-2. (6) 

For ions entering the sheath with the Bohm velocity mq = — 1 . The field at the effec- 
tive wall as a function of the wall potential Tj„, — rj ((^o) is then given by 



I7vv - - \/2V2tJvv + 1 + 2exp(-77„) - 4. 



(7) 
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The total surface density of ions in the sheath equals the total surface density of 
electrons in the ESL A^. It can be related by Poisson's equation to the field at the 
wall and reads 

^ - = = ^^(zo) = -noXo,:.. (8) 

Combing Eqs. ([8]) and (j7]l gives the total surface density of electrons to be inserted 
into the ESL as a function of the wall potential. 

The wall potential itself is determined by the flux balance condition, = j,, 
which, in the ESL model, is assumed to be fulfilled at z = zq. Using the Bohm flux 
for the ions and the thermal flux for the electrons, 

the wall potential is given by HI 
that is, 

^.^-^^Inf^V (11) 

In the collision-less sheath model the wall potential depends only on the electron 
temperature and the ion to electron mass ratio. 
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Fig. 5 Effective potential for the graded interface on which the model of an ESL is based. It 
encompasses the image potential as well as the potential offset at the interface. 



We can now move on to the distribution of surplus electrons in the surface layer. 
It is primarily determined by the potential in the interface region surrounding an 
ideal dielectric interface. We first single out the relevant microscopic potentials at 
the interface (illustrated in Fig. |4]l and then describe an effective model for them 
(Fig.ig. 

Far away from the surface the potential is given by the repulsive sheath potential. 
Closer to the surface an attractive surface potential, the image potential, which stems 
from the dielectric mismatch at the interface, becomes dominant. Classically the 
image potential is given by the expression |25 1 



0im(z) 



4(e.v + l)i 



(12) 



where is the static dielectric constant. This expression is valid far away from the 
wall. Closer to the wall the 1 /z shape is modified as the potential merges continu- 
ously with the crystal potential in the bulk of the dielectric. Even more important 
inside the dielectric is the band structure. For a dielectric the valence band is full 
while the conduction band is completely empty. Thus surplus electrons can popu- 
late the conduction band and the conduction band minimum acts as a long range 
potential for the electrons. In addition to the image potential, the surface potential 
comprises an offset due to a charge double layer formed immediately at the surface. 
It is due to electron density leaking out of the solid as a consequence of free energy 
minimisation or to relaxation and polarisation of the atomic bonds at the surface as 
a consequence of the truncation of the lattice. This offset of the conduction band 
minimum at the surface to the potential in front of it is characterised by the elec- 
tron affinity of the solid %. For % <Q the conduction band minimum lies above the 
potential just in front of the solid, for j > it lies below it. So far we have consid- 
ered an uncharged surface. This raises the question whether the electronic structure 
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of the surface is changed by the surplus electrons from the plasma. In comparison 
to the electrons responsible for the chemical binding within the dielectric the ad- 
ditional electrons coming from the plasma are only a few. The available electronic 
states and the offset of the energy bands in the bulk with respect to the potential out- 
side the dielectric will thus not be changed significantly by the presence of the wall 
charge. Note however that a chemical modification of the plasma wall can change 
the electronic structure. In particular the electron affinity is susceptible to surface 
coating and adatoms. Surface termination by elements with small electronegativity, 
for instance hydrogen, may induce a negative electron affinity 1 26 1 while elements 
with larger electronegativity, for instance oxygen, can lead to a positive electron 
affinity ll27l . 

From these microscopic potentials we construct an effective surface potential. 
To obtain a realistic image and offset potential without performing an atomisti- 
cally accurate calculation we employ the model of a graded interface. This model 
is parametrised by experimentally measured (or where not available theoretically 
calculated) values for the electron affinity, the dielectric constant, and the conduc- 
tion band effective mass. First proposed by Stern 1211 to remove the unphysical 
singularity of the image potential at z = and later extended to potential offsets at 
semiconductor heterojunctions 1281 it assumes the smooth variation of parameters 
that change abruptly at the interface over a length on the order of the lattice constant. 
For our purposes we assume a sinusoidal interpolation of the dielectric constant, off- 
set of the long range potential and effective electron mass with a grading half-width 
of d ^ SA |20|. The model does not account for effects associated with intrinsic 
surface states (Shockley or Tamm states |29|) or additional states which may arise 
from the short-range surface potential. Nevertheless the graded interface model is a 
reasonable description of the surface, suited for dielectrics with ionic bonds which 
typically have no intrinsic surface states. 

The total surface potential 

0surf(z) = ^im{z) + 0offset(z) (13) 

comprises the graded image and offset potential. It is continuous across the crys- 
tallographic interface at z = and enables us thereby to also calculate a smoothly 
varying electron distribution in the ESL. 

Using Eq. ([TJ we can now determine the position zo of the effective wall, that is, 
the maximum extent of the ESL on the plasma side. The derivative of the surface 
potential is 0^'^^.^ — ^og^gj + Due to the relatively weak field in the sheath com- 
pared to the image force, the boundary zo will be so far away from the interface that 



'^'offset vanishes and the image potential obeys ( 12 1. Thus, the boundary between the 
ESL and the plasma sheath is given by 



with (pl. = -{kBTerii,)/{eXD), and 7]^. given by 
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We can now turn to the distribution of the plasma-supplied excess electrons in 
the ESL. For this quasi-static model we assume that the wall-bound surplus elec- 
trons are in thermal equilibrium with the wall. Thus, their distribution minimises 
the grand canonical potential and satisfies Poisson's equation. Inspired by Tkharev 
and Danilyuk Ii22l we apply density functional theory ll30l |3TI to the graded inter- 
face. For the purpose of this exploratory calculation we will use density functional 
theory in the local approximation. Quite generally, the grand canonical potential of 
an electron system in an external potential V (r) is given in the local approximation 
by 



Q.^ y y(r)«(r)dr- -y 0c(r)«(r)dr + G[«]-juy «(r)dr, (15) 

where G[n] is the grand canonical potential of the homogeneous system with density 
n{r) and the Coulomb potential is determined by 

V (£(r)V0c(r)) = ATten{r). (16) 

The ground state electron density minimises £2, that is, it satisfies 

V{r)-e<^c{r)+ll\n)-^=Q, (17) 

where /x''(n) — 5G[n\/5n is the chemical potential for the homogeneous system. 



In our one-dimensional graded interface model Eq. ( 17 1 reduces to 

-e<^[z)+^\z)-^=0, (18) 
where /i''(z) = jJ.'^{n{z),T), and the electrostatic potential 

^{z) = (j>.uAz) + ^c{z) (19) 

consists of the potential of the bare surface given by Eq. ( pj) and the internal 
Coulomb potential which satisfies Poisson's equation 

^(e{z)^^c{z)) ^Anen{z) (20) 



dz \ dz 

with the graded dielectric constant e(z). The boundary conditions 0c (zo) = and 
(l>c{zo) = guarantee continuity of the potential at zo and include the restoring 
force from the positive charge in the sheath. Note that the Coulomb potential derived 
from this equation already includes the attraction of an electron to the image of 
the charge distribution. The image potential contains only the self-interaction of an 
electron to its own image. 

For the functional relation /i''(z) = /i''(n(z), T) we take the expression adequate 
for a homogeneous, non-interacting, non-degenerate electron gas 
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Fig. 6 Plasma-supplied excess electron density n (upper panel) and the potential — ^ (lower panel) 
it gives rise to for an LiF (left) and AI2O3 (right) surface in contact with a helium discharge with 
fio = 10^ cm^-' and k^T^ = 2eV. 



where m(z) is the graded effective electron mass. This is justified because the density 
of the excess electrons is rather low and the temperature of the surface is rather high, 
typically a few hundred Kelvins. 

For the calculation of the quasi-stationary distribution of surplus electrons Eqs. ( fTSj l 
and ( pOl have to be solved iteratively (until is stationary) subject to the additional 
constraint j^" dz«(z) — N which guarantees charge neutrality between the ESL and 
the plasma sheath. The position Zs < is a cut-off which has to be chosen large 
enough in order not to affect the numerical results. In a more refined model for the 
ESL the crossover of the wall charge to the neutral bulk of the dielectric can be 
taken into account by splitting the ESL into an interface specific region and a space 
charge region. This approach which turns Zs from an ad-hoc cutt-off into the bound- 
ary between the two regions is described in |20|. For the charge distribution at the 
interface the simple model is however sufficient. 

We now use the ESL model to calculate for a helium discharge in contact with 
a LiF and AI2O3 surface the electron density and potential across the plasma wall. 
Our main focus lies in the identification of generic types of electron distributions in 
the ESL. 

Depending on the electron affinity one of two scenarios for the distribution of 
surplus electrons in the ESL may be realised. For negative electron affintiy J < 0, 
the conduction band minimum lies above the long-range potential just in front of 
the surface. It is thus energetically unfavourable for surplus electrons to populate 
the conduction band. Instead they are located in the image potential in front of the 
surface. The energy of an electron in the image potential reaches its minimum just in 
front of the surface at the beginning of the graded interface. This is confirmed by the 
electron density and the potential for LiF (x < 0) shown in the left panel of Fig. |6] 
For X < 0, the excess electrons coming from the plasma form therefore an external 
surface charge in the image potential in front of the crystallographic interface. As 
the electrons do not penetrate into the soUd the band bending associated with it is 
negligible. The external surface charge is very narrow. It can thus be considered as 
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a quasi-two-dimensional electron gas, similar to the surface plasma anticipated by 
Emeleus and Coulter ifTSllSZI. 

For positive electron affinity X > ^ the situation is dramatically different. In this 
case the conduction band minimum is below the potential just outside. It is thus 
energetically favourable for electrons to accumulate inside the dielectric. This can 
be seen in the right panel of Fig.|6]which shows the electron density and the potential 
for an AI2O3 surface (x > 0). Although the image potential still induces an attractive 
well in front of the surface the minimum of the potential energy for electrons — e0 
is reached inside the dielectric. Excess electrons coming from the plasma are thus 
mostly located inside the wall where they form a space charge in the conduction 
band. The electron distribution extends deep into the bulk and entails bending of the 
energy bands near the surface. Note the different scales of the axes for the left and 
right panels of Fig.|6] On the scale where variations in the space charge for AI2O3 
are noticeable the electron distribution for LiF is basically a vertical line. 

So far we have shown the potential and the electron distribution in the ESL. Now, 
we will compare the potential and the charge distribution in the ESL with the ones 
in the plasma sheath. The wall bound electrons in the surface layer constitute an 
electron distribution in thermal equilibrium with the wall which balances electron 
and ion influx at the sheath-ELS boundary zo- On top of this quasistationary elec- 
tron distribution there are fluxes of electrons and ions, coming from the plasma. 
These currents, which begin at zo and persist to the location in the ESL where ei- 
ther electron trapping or electron-ion recombination takes place are not described 
by our simple ESL model. The electron and ion densities in this model are thus 
discontinuous at zo- The potential, however, which has been obtained from the inte- 
gration of Poisson's equation is continuous and differentiable everywhere. Between 
the crystallographic interface and zo the electron and ion flux from the plasma would 
be important. The neglect of the charge densities associated with these fluxes does 
however not affect the potential because they are too small to have a significant 
effect. 

Figure|7]compares the ESL and the plasma sheath for LiF (x < 0, left) and AI2O3 
(X > 0, right). We plot the electron and ion densities (upper panel) as well as the 
electric potential (lower panel). Note the logarithmic scale in the left panel and 
the linear scale in the right panel. For LiF the surface electrons are bound in the 
image states in front of the surface due to the negative electron affinity. Far from 
the surface, the potential approaches the bulk plasma value chosen to be zero. In the 
sheath the potential develops a Coulomb barrier and reaches the wall potential (j)^. 
at Zo, the distance where the sheath merges with the surface layer (vertical dotted 
line). The wall potential is the potential just outside to which the energies of the bulk 
states are referenced. Closer to the surface the potential follows the attractive image 
potential while at the surface the repulsive potential due to the negative electron 
affinity prevents the electron from entering the dielectric (only scarcely seen on the 
scale of the figure). 

For AI2O3 the excess electrons constituting the wall charge penetrate deep into 
the dielectric due to the positive electron affinity. Compared to the variation of the 
electric potential in the sheath the band bending in the dielectric induced by the 
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Fig. 7 Density of plasma-supplied surplus electrons trapped in the ESL, electron and ion density 
in the plasma sheath, and electric potential for a LiF (left) and an AI2O3 surface (right) in contact 
with a helium discharge with plasma parameters no = lO'cm^' and ksT^ = 2eV (taken from |33 1). 
The crystallographic interface is at z = 0. Note the different scales of the two panels. The deep 
penetration of the AI2O3 wall charge is due to the neglect of defect states and other collision 
centres. 

wall charge is rather small as indicated by the variation of inside the dielectric. 
This is because is large and most surplus electrons are bound close to the crys- 
tallographic interface. The interface specific region where the short ranged surface 
potential varies strongly is very narrow compared to the extent of the two space 
charge regions. 

The ESL should be regarded as the ultimate boundary layer of a bounded gas 
discharge. Our model for the ESL captures the coupling of the neutral bulk plasma 
to an electron plasma at the boundary separated from the bulk plasma by the plasma 
sheath. Depending on the electron affinity this boundary plasma can be either a 
quasi-two-dimensional electron plasma or an electron(-hole) plasma inside the wall. 



3 Electron physisorption 

We can now turn to the electron kinetics in the ESL. Depending on the electron 
affinity two distinct scenarios for electron trapping emerge which are schematically 
shown in Figs. [8] and [9] For negative electron affinity (see Fig. [8]l electrons that 
overcome the sheath potential do not penetrate the solid as their energy falls in the 
band gap where no internal states are available. Instead, they are trapped in the im- 
age potential in front of the surface. Energy-relaxation is in this case due to surface 
vibrations which trigger transitions to the bound image states. As the image poten- 
tial is relatively deep compared to the energies of the longitudinal acoustic phonon 
responsible for electron trapping and the wavefunctions of the unbound states are 
suppressed close to the surface the probability for sticking is small. This scenario, 
which we have investigated in detail Il34l[35ll36]| . will be presented below. 
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Fig. 8 Electron trapping in the ESL in front of a dielectric with % 



For positive electron affinity (see Fig. [9]l electrons that overcome the sheath po- 
tential penetrate into the solid. There they initially occupy high-lying states of the 
conduction band. Subsequent electron energy relaxation is due to scattering with 
bulk phonons of the dielectric and the plasma-supplied electrons are eventually 
trapped at the bottom of the conduction band in the potential well created by the 
restoring force from the plasma. Due to the high-density of states and an efficient 
scattering mechanism we expect the sticking coefficient to be much larger in this 
case. However, a detailed analysis of this scenario still remains to be done. 

We will now consider the physisorption of an electron in the image potential in 
front of a dielectric with < in detail. We describe the time evolution of the occu- 
pancy of the bound surface states with a quantum-kinetic rate equation fjT', '381. It 
captures all three characteristic stages of physisorption: initial trapping, subsequent 
relaxation and desorption ||39ll401 . 

The time dependence of the occupancies of the bound states is given by Il37ll38l 

-^n„(f) =£[W„„/n„'(0-^«'«n«(0]-L^ten«(0+LW«U-l(f) (22) 
n' k k 

which can be rewritten as 

4nn(f) =Y.Tnn'nn'{t)+Y,^tWnkjk{t), (23) 

n' k 

where W„i„ is the probability per unit time for a transition from a bound state n to 
another bound state n', Wii„ and W„k are the probabilities per unit time for a transition 
from the bound state n to the continuum state k and vice versa, and Tf = 2L/v- is 
the transit time through the surface potential of width L, which, in the limit °o, 
can be absorbed into the transition probability. The matrix r„„/ is defined implic- 



itly by the above equation. The last term in Eqs. (22 1 and (23 1, respectively, gives 
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Fig. 9 Electron trapping in the ESL in a dielectric with X > 0. 



the increase in the bound state occupancy due to the coupHng between bound and 
unbound surface states. 

The probability for an approaching electron in the continuum state k to make a 
transition to any of the bound states is given by the prompt energy-resolved sticking 
coefficient 

^tj^^nk- (24) 



prompt _ 
^e,k — 



Treating the incident electron flux as an externally specified parameter, the so- 



lution of Eq. ( 22 1 describes the subsequent relaxation and desorption. It is given 
by 



(25) 



where e„ and e„ are the right and left eigenvectors to the eigenvalue — A^- of the 
matrix T„„i. 

If the modulus of one eigenvalue, Ao, is considerably smaller than the moduli 
of the other eigenvalues, A^-, a unique desorption time and a unique sticking coef- 
ficient can be identified |38|. In this case Ao governs the long time behaviour of 
the equilibrium occupation of the bound states, «^'' ^ g-EnlksTs^ inverse can 

be identified with the desorption time, Aq"' = Te- The bound state occupancy «„(f) 
splits into a slowly varying part n^{t) given by the K — Q summand in Eq. (25 i and 
a quickly varying part n{ (f ) given by the sum over (f 7^ in Eq. (25 1. 



The electron adsorbate, i.e. the fraction of the trapped electron remaining in the 
surface states for times on the order of the desorption time, is given by the slowly 
varying part only, n^[t) = Differentiating «°(f) with respect to the time. 



dt 



n\t)^l^s^r^jk{t)-:k^n\t) 



(26) 
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unbound states 



Fig. 10 Schematic represen- 
tation of the series of image 
states in the image potential. 
The lowest bound states lies 
very deep and is connected 
to a group of upper bound 
states only by a two or three- 
phonon processes. The upper 
bound states are connected by 
one-phonon processes to the 
continuum. 



2-3 phonon deep 




upper bound states 



lowest bound state 



we can, following Brenig BOl . identify the kinetic energy-resolved sticking coeffi- 
cient 

*" = ^rL4''^~"°^^«^ (27) 

n.n' 

giving the probability for both, initial trapping and subsequent relaxation. 

If the incident unit electron flux corresponds to an electron with Boltzmann dis- 
tributed kinetic energies, the prompt or kinetic energy-averaged sticking coefficient 
is given by 

^-W^' '''' 

where jS^' = kgTe is the mean electron energy. 

These kinetic equations rely on the knowledge of the transition probabilities. 
They have to be calculated from a microscopic model for the electron-surface in- 
teraction. The image potential supports a series of bound states which can trap 
electrons temporarily at the surface. Transitions between image states are due to 
dynamic perturbations of the surface potential. The image potential is very steep 
near the surface. A particularly strong perturbation arises therefore from the surface 
vibrations induced by the longitudinal acoustic bulk phonon perpendicular to the 
surface. Describing them with a Debye model, the maximum phonon energy of this 
mode is the Debye energy h(OD- Measuring energies in units of the Debye energy, 
important dimensionless parameters characterising the potential depth are 

F F —F I 

£„ = — — and A„„, = —— — , (29) 
ncoo n(£)o 

where £„ < is the energy of the n* bound state. 

The depth of the image potential can be classified with respect to the Debye 
energy. If + 1 > Zii2 > — «, we call the potential «-phonon deep. If ei > — 1 we 
call it shallow. For a shallow potential the lowest bound state is coupled by a one- 
phonon transition to the continuum. If the potential is n-phonon deep an n-phonon 
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process links the lowest bound state to the second bound states, which is linked by 
one-phonon transitions either directly or via intermediates to the continuum. For the 
dielectrics which we consider the potential is two-phonon (MgO) or three-phonon 



deep (LiF). A schematic representation of the image states is given in Fig. 10 

For the calculation of the transition probabilities used in the rate equation we 
need a microscopic model for the image potential and the electron-surface vibration 
interaction. Classically the image potential ensues from the dielectric mismatch at 
the interface. On a microscopic level it arises from the coupling of the electron to a 
polarisable surface mode of the solid. For a dielectric this is a surface phonon. For 
LiF and MgO the low frequency dielectric function is dominated by a transverse 
optical (TO) phonon with frequency 0)7-0. It can be approximated by BTl 

e(to) = l + (e,-l)^^^, (30) 

where is the static dielectric constant. The bulk TO-phonon gives rise to a sur- 
face phonon. Its frequency = Wro + £s) /2 is determined by the condition 
e{(Os) — — I. The electron couples to this surface phonon according to Ii41ii42j 



H-- — 
2m 



A+Y, fiO)sal aq + £ (^q (x,, )flQ + (pX (x,, ) (31) 
Qii Qii 



with 



where Aq — {Es — 1)/ (4(ej + 1)) and the subscript || denotes vectors parallel to the 
surface and Oq^J are annihilation (creation) operators for surface phonons. 

Applying a unitary transformation |42| separates the coupling into a static part 
which takes the classical form of the image potential ^ l/z and a dynamic part of the 
electron-surface phonon coupling, which encodes recoil effects and encompasses 
momentum relaxation parallel to the surface. While the classical image potential al- 
lows a simple description of the image states which captures their properties fairly 
well it is not sufficient for the calculation of probabilities for surface vibration in- 
duced transitions, as they perturb the electronic states most strongly close to the sur- 
face where the image potential is steepest. Unfortunately, in this region the classical 
image potential has an unphysical divergence. To remedy this we use a variational 
procedure to extract the static image potential from 1411. Thereby we keep some 
recoil effects which make the recoil-corrected image potential ^ l/{z + Zc) with a 
cut-off parameter Zc = ^JhjTmO)^ |34| divergence free. 

Transitions between the eigenstates of the recoil-corrected image potential are 
due to the longitudinal-acoustic bulk phonon perpendicular to the surface. The 
Hamiltonian from which we calculate the transition probabilities is given by 



(33) 
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where 

i/f'-^^V^c, (34) 

describes the electron in the recoil-corrected image potential, 

Hph=Y,rzcoQblbQ (35) 

e 

describes the free dynamics of the bulk longitudinal acoustic phonon responsible for 
transitions between surface states, and 



H'i;,^j:{qKiu,z)\q)cl,c, (36) 



denotes the dynamic coupling of the electron to the bulk phonon. 

The perturbation Vp{u,z) can be identified as the difference between the dis- 
placed surface potential and the static surface potential. It reads, after the transfor- 
mation Z — Zc, 

Vp[u,z) = ■ \ . (37) 

z + u z 

In general, multi-phonon processes can arise both from the nonlinearity of the 
electron-phonon coupling H'^^^i^ as well as from the successive actions of //gl^/j 
encoded in the T matrix equation, 

T-H^'-p„+H^'-UGoT, (38) 



where Go is given by 



Go= (£-//f"=-//p/, + /0+) (39) 



The transition probability per unit time from an electronic state q to an electronic 
state q' encompassing both types of processes is given by 1.43 J 

W,>,= —Y.^^-^\{s',q'\T\s,q)\^5{E,-E,+E,-E^,), (40) 



s,s' 



where j3j = {kgTs)^^, with the surface temperature and \s) and \s') the initial 
and final phonon states. We are only interested in the transitions between electronic 
states. It is thus natural to average in Eq. (40 1 over all phonon states. The delta 
function guarantees energy conservation. 

In principle, multiphonon transition rates can be obtained by iterating the T- 
matrix and evaluating Eq. (40l. Up to ^(m^), for instance, the T-matrix reads 



T = V1+V2 + V3+V1G0V1 +V2GQV1 + V1G0V2 + V1G0V1G0VU (41) 
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Fig. 11 Left: Inverse desorption time T^^' in the two-phonon approximation fox Tq/T = 5 as a 
function of the Debye temperature To- Data for To = 2500^ apply to graphite. Right: Prompt 
energy-resolved sticking coefficient as a function of the electron energy for a two phonon deep po- 
tential (To = 2500Ar and = "iSlK). The contribution of one-phonon processes (red) far outweighs 
the contribution from two phonon processes (blue). 



where the Vi ^ u' originate from expanding Eq. (37 1 in the displacement u. The T- 
matrix enters as {T){T*) into the transition probabiHty. The term (Vi)(V'j*) can be 
identified as the Golden Rule transition probability. Proportional to it is a one- 
phonon process. Two-phonon processes, proportional to m^, are represented by the 
terms 

{Vi){V,*),{V3){V:)AVi){V^GlV,*), {V,GoV2){V,*),{Vi){V,*GlV^), 

{V2GoVi){V{),{Vi){V{GlV*GlV*),{ViGoViGoVi){Vl), (42) 



and 



(Vi) (V2*> , (ViGoVi) (v;G5yr) , {V2) (vrCSv;) , {ViGoVi){V^). (43) 



A complete two-phonon calculation would take all these processes into account as 
they stand. This is however not always necessary. A closer analysis (see 1341 ) of 
the first group of terms reveals that they are two-phonon corrections to transitions 
already enabled by a one-phonon process. We assume that theses corrections are 
small and evaluate only the second group for transitions where they enable two- 
phonon transitions which are not merely corrections to one-phonon transitions. The 
details of the evaluation of the transition probabilities, including a regularisation of 
divergences by taking a finite phonon lifetime into account, can be found in If34ll35]| . 

The expansion of the T-matrix allows the calculation of transition probabilities 
for the two-phonon deep potentials of graphite and MgO. However, for a three- 
phonon deep potential, for instance for LiF and CaO, this approach is no longer 
feasible. From Refs. Il34ll35]| we qualitatively know the relevance of different types 
of multi-phonon processes. For continuum to bound state transitions, for instance, 
one phonon processes are sufficient at low electron energies. We will therefore com- 
pute the transition probabilities between bound and continuum states in the one- 
phonon approximation. For transitions between bound states, we found that multi- 
phonon processes due to the nonlinearity of the electron-phonon couphng tend to 
be more important than the multi-phonon processes due to the iteration of the T- 
matrix. Hence we expect that an approximation which takes only the nonlinearity of 
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Fig. 12 Energy resolved 
prompt sticking coefficient 
for MgO (left) and schematic 
representation of the bound 
and continuum wavefunctions 
(right). The suppression of the 
amplitude of the continuum 
wavefunction close to the 
surface for small electron 
energies causes the small 
sticking coefficient. 




n=3 E=-45meV 



the electron-phonon interaction nonperturbatively into account to be sufficient for 
the identification of the generic behaviour of muhi-phonon-mediated adsorption and 
desorption. This approach is described in ll36l . 

We now turn to resuhs for the desorption time and the sticking coefficient for 
physisorption in the image states. Before showing results for MgO and LiF, we con- 
sider the dependence of the electron kinetics on the potential depth and the relevance 
of one or two-phonon processes. For this we show in the left panel of Fig. [TT|the 
inverse of the desorption time as a function of the Debye temperature To — hcoo/kB 
which sets the energy scale of the acoustic phonons. While the absolute depth of the 
potential remains constant, varying the phonon energy tunes the effective potential 
depth. Fig.[TT|reveals that for a shallow potential desorption is most efficient while 
for a one-phonon deep potential it gradually becomes less efficient. When the po- 
tential becomes two-phonon deep desorption suddenly becomes even slower which 
reflects the small magnitude of two-phonon transitions compared to one-phonon 
transitions. This justifies our approximation of neglecting two-phonon corrections 
to transitions already enabled by a one-phonon process. For a Debye temperature 
Tj) — 2500^r the results apply to graphite. 

The energy-resolved prompt sticking coefficient, plotted in the right panel of Fig. 
1 1 for graphite, shows that also for bound state-continuum transitions one-phonon 
processes far outweigh two-phonon processes. Thus, sticking is mainly due to one- 
phonon processes. Moreover we find that the sticking coefficient drops sharply at 
specific energies. These accessibility thresholds occur whenever one bound state 
becomes no longer accessible from the continuum because the energy difference 
exceeds one Debye energy. 

The sticking coefficient takes relatively small values — on the order of 10^^. 
This is due to the small matrix elements between bound and continuum states. 
Fig. [12] shows the wave functions of representative bound and continuum states 
for MgO. For large energy the continuum wave functions have sinusoidal shape 
whereas for small energies their amplitude is significantly suppressed close to the 
surface. These two behaviours correspond to two limits for the wavefunction. For 
simplicity we discuss this for a 1 /x potential without cut-off. The continuum wave- 
functions read ^ki^) ^ ~7t^-i ' where x = z/ob- For x — > oo we ob- 
tain (pk{x) ^ sin(Aofcc) which also holds for large k while for A; we obtain 
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Fig. 13 Inverse desorption time Tj, ' (left), prompt energy resolved sticking coefficient (middle) 
and prompt as as well as kinetic energy averaged sticking coefficient (right) for MgO and LiF. 



^k{^) ^ Vky/2AoxJi {y/SAox). The proportionality to \/k entails a strong suppres- 
sion of the wavefunction for low energy. 

Electron physisorption at a dielectric surface with negative electron affinity is an 
intriguing phenomenon due to the interplay of potential depth, magnitude of ma- 
trix elements and surface temperature 1341 [35l U6 1 . Initial trapping of the electrons, 
characterised by the prompt sticking coefficient, occurs in the upper bound states 
by one-phonon transitions. Relaxation after initial trapping depends on the strength 
of transitions from the upper bound states to the lowest bound states. If the lowest 
bound state was linked to the second bound state by a one-phonon process the elec- 
tron would relax for all surface temperatures. If they are linked by a multiphonon 
process relaxation takes place only for low temperature, whereas for room temper- 
ature a relaxation bottleneck ensues as the electron desorbs from the upper bound 
states before it drops to the lowest bound state. This leads to a severe reduction of the 
kinetic sticking coefficient compared to the prompt sticking coefficient. The dom- 
inant desorption channel depends also on the depth of the potential. For a shallow 
potential desorption occurs directly from the lowest bound state to the continuum. 
For deeper potentials desorption proceeds via the upper bound states. Desorption 
occurs then via a cascade in systems without and as a one-way process in systems 
with relaxation bottleneck |36|. An overview of results for the desorption time as 
well as prompt and kinetic sticking coefficients for MgO and LiF is given by Fig. 
13 Most important for the plasma context is that Se < 1 and T^' ^0 Ii34l[35ll36ll . 
implying that a dielectric surface with a negative electron affinity is not a perfect 
absorber for electrons. 



4 Mie scattering by a charged particle 

In the previous section we have employed our model for the ESL to analyse the elec- 
tron kinetics at the plasma boundary and to calculate desorption times and sticking 
coefficients for a dielectric with negative electron affinity. Another application of the 
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model for the ESL is the study of the optical properties of the wall-bound electrons. 
In this section we will study the effect of surplus electrons on the scattering of light 
by a sphere (Mie scattering) with an eye on a charge measurement for dust particles. 
Most established methods for measuring the particle charge 1441 1451 |46l require 
plasma parameters which are not precisely known. The particle charge, however, 
is the fundamental quantity of interest in the field of dusty plasmas. It determines 
the interactions between dust grains and between dust grains and external electro- 
magnetic fields. While Mie scattering is akeady used in the field of dusty plasmas 
to determine the particle size |47 |, it has up to now not been used to determine the 
particle charge. Quite generally the Mie signal contains information about excess 
electrons as their electrical conductivity modifies either the boundary condition for 
electromagnetic fields or the polarizability of the material ll48l l49l l50l [5 1 J . But to 
what extent and in what spectral range the particle charge influences the scattering 
of light remains an open question. 

The location of the surplus electrons depends on the electron affinity % of the 
dust grain. For < 0, as it is the case for MgO, CaO or LiF Il52l l53 l, electrons are 
trapped in the image potential in front of the surface where they form a spherical 
two-dimensional electron gas around the grain. The image potential results from the 
static interaction of the electron to a surface mode associated with the transverse 
optical (TO) phonon. The residual dynamic interaction limits the surface conduc- 
tivity Os of the image bound electrons. The electrical field parallel to the surface 
induces a surface current K = a,E|| which is proportional to the surface conduc- 
tivity Os- The surface current enters the boundary condition for the magnetic field. 
Thus, X (H; + - H,) = ^K, where c is the speed of light 

For > 0, as it is the case for AI2O3 or Si02, electrons form a space charge 
layer in the conduction band 1201 . Its width, limited by the screening in the dielec- 
tric, is typically larger than a micron. For micron sized particles we can thus assume 
a homogeneous distribution of the excess electrons in the bulk. The effect on light 
scattering is now encoded in the bulk conductivity of the excess electrons Of,, which 
is limited by scattering with a longitudinal optical (LO) bulk phonon. The bulk con- 
ductivity gives rise to a polarizability a = 'XniOi, j CO, with CO the frequency, which 
enters the refractive index = ^/e + a. 

To solve the scattering problem we expand the incident plane wave, 

E,- = EQe''''-""e„ H, = -^oe'^-'^'^'ev (44) 

CO 

in spherical vector harmonics memn,^omn, T^emn, and no,„„ |!54l which are solution of 
the vector wave equation V^C + fc^C=0 , in spherical coordinates 

E-=^o.-'-E."j^K!,-,„i;!V (45) 



,S n(n + l) 



ck 



H,- = --E,e-'-' £ (mill +/nW ) (46) 

CO «(n + l) V 
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and match the incident wave with the reflected wave, 

E, = Eoe-'-' £ i"^^^ (<mf )„ - /^>i^i) , (47) 



n(n + l) 



and the transmitted wave. 



rt=l V ' / 



n=l 



(^ is the wavenumber inside the particle) at the boundary of a dielectric sphere with 
radius a. The boundary conditions at the surface are given by 

CrX (C, + C,-C,) =0 (51) 

for C = E, H in the case of j > and in the case of j < by the above equation for 
E and 

An 

e, X (H, + H, - H,) = — (T.E|| . (52) 

c " 

This gives, following Bohren and Hunt IM |, the scattering coefficients 

r ^ Vn {Np )V„{P)- [NW'„ {Np ) - /T v/„ (Np)] y/,, (p ) 

" [NxifUNp)-nxif„{Np)]Up)-Vn{Np)^l,{p) ^ ' 

and 

.r ^ ¥n (Np ) y/,, (P ) - [NXtfn (Np ) + /T < (Np)] Xtf'„ (p ) 

« [iVv^„(A^p) + /TV/UA^p)]^^(p)-r;;(A^p)^«(p) ^ ' 

where for J < (j > 0) the dimensionless surface conductivity t{(o) ~ 4nGs{co)/c 
(t = 0) and the refractive index — ^/e (N = ^/E~+a), p = ka = Inajk the 
size parameter , where k is the wavenumber, y/nip) — \J ^p /^Jn+i/iiP)^ ^«(p) — 
yjnp l2H^^^^ir^{p) with Jn{p) the Bessel and Hj,^\p) the Hankel function of the 
first kind. As for uncharged particles the scattering efficiency is given by 

Qs = ^i{2n+l){\a';,\HK\^), (55) 

P n=l 

the extinction efficiency by 

Qr = -^t{2n + l)3i{a:+b:) (56) 
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and the absorption efficiency by Qa — Qt ~ Qs- Any effect of the surplus electrons 
on the scattering of light, encoded in and b"^, is due to the surface conductivity 
(X < 0) or the bulk conductivity {% > 0). 

For X <Owe use a planar model, which we justify below, to calculate the surface 
conductivity. For the calculation of the surface conductivity we will employ the 1 /z 
image potential. The electron-phonon interaction is given by Eq. pT| ). We apply 
the unitary transformation H — > UHU^^ with U — e'^, S — {i/fi(Os)Y.Q^^ '^Q\\^Q\\ ~ 

oL <1>Q to separate static and dynamic coupling B2ll . The former leads to the image 
^11 ^11 

potential V ~ — Aoe^ / z with Aq = (eo — 1 ) / (4 (eo + 1 ) ) supporting a series of bound 
Rydberg states whose wave functions read 



^«u(x,z)^^.-J^W„,/J?^) (57) 
VA V aBnnl-^ ' \ nag ) 

with flfi the Bohr radius, k = (kx^ky), x = (x,y), A the surface area and Wn.m Whit- 
taker's function |55|. Since trapped electrons are thermalised with the surface and 
the spacing between Rydberg states is large compared to k^T , they occupy almost 
exclusively the lowest image band n — I. Assuming a planar surface is justified 
provided the de Broglie wavelength X^ig of the electron on the surface is smaller 
than the radius a of the sphere. For a surface electron with energy E\r^a/kB = SOO/T 
one finds X^b ~ 8 x 10^^ cm. Thus, for particle radii a > 10 nm the plane-surface 
approximation is justified. The residual dynamic interaction enables momentum re- 
laxation parallel to the surface and hence limits the surface conductivity. Introduc- 
ing annihilation (creation) operators c^^' for electrons in the lowest image band, the 
Hamiltonian describing the dynamic electron-phonon coupling in the lowest image 
band reads ll56l 

H = Y^ ekc^ck + h 0), flgaQ + Hin, (58) 

k Q 

with ^ 

Hint = ^ X!^k,Q4+Q («Q^«-q) (^9) 

vAkQ ^ ' 

where the matrix element, calculated with the wave function given by Eq. pT] ), is 



MkQ = 



le^/iiA^ ( 2Ao 



eflB + 2Ao 



(60) 



(m is the electron mass). Within the memory function approach f57l the surface 
conductivity can be written as 
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_ t' (OK) — t' (300K) . . . t" (free e) — -a' (OK) — -a' (300K) ... -a' (free e) 




Fig. 14 Dimensionless surface conductivity T = t' + it" for LiF for = 10'^ cm^^ (left) and 
polarisability of excess electrons a = a' + ia" for AI2O3 for = 3 x lO'^ cm^' (right) as a 
function of the inverse wavelength A ^ ' . 



with «s the suiface electron density. The memory function is then evaluated up to 
second order in the electron phonon coupling ll58l . Since M(a)) is independent of 
« v the surface conductivity is proportional to the surface density of electrons n, . 

For X > the interaction of the electron with a longitudinal optical bulk phonon 
limits the bulk conductivity. The coupling of the electron to this mode with fre- 
quency (Olo is described by ||59l 

= ^ L T'^k+q^k (flq + fl-q) , (62) 
V V k,q ? ^ ^ 



where M — Ine'^TKOw (£» ' ^ £0 0- ^'^^ calculation of the bulk conductivity 
we employ again the memory function approach. In this case the bulk conductivity 
is proportional to the electron density n^. 

The charge effect on scattering is encoded in the surface conductivity for % <^ 
and the bulk conductivity for ;|; > which enter through the dimensionless sur- 
face conductivity T or the polarisability a into the scattering coefficients. Figure [l4] 
shows 1 — x' + ix" for LiF and the polarisability for a — a' + a" for AI2O3 as a 
function of the inverse wavelength A ^ ' . They turn out to be small even for a highly 
charged particle with = lO'^ cm^^ (corresponding to 3 x lO'^ cm^^ for 
X > and a = 1/im). Compared to a free electron gas where M(a)) = (implying 
t' = and a" = 0), the electron-phonon coupling reduces t" and a' considerably. 
For T = OK,t' =0 {a" = 0) for A"' < A,r' = 675 cm"', the inverse wavelength of 
the surface phonon ( A"' < = 807 cm"', the inverse wavelength of the bulk LO 
phonon) since light absorption is only possible above the surface (bulk LO) phonon 
frequency. At room temperature t" and a' still outweigh t' and a". The tempera- 
ture effect on z" is less apparent for A"' > 300 cm"' than for a' but for A"' < 300 
cm"' a higher temperature lowers t" considerably. 
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Fig. 15 Dielectric constant e = e' + /e", refractive index N = n + ik (first row) and extinction 
efficiency Q, (second row), scattering efficiency Qs (third row) as well as absorption efficiency Qa 
(fourth row) depending on the particle radius a for LiF, CU2O and AI2O3 as a function of A^' . 



We now turn to the scattering properties of the sphere. For LiF and AI2O3 one 
or two phonon modes dominate the dielectric constant. Far above the highest TO 
phonon mode (A^' > 700 cm^ ' for LiF and A"' > 900 cm for AI2O3) the re- 
fractive index is real, relatively small and does not vary strongly with frequency. 
In this regime a micron sized grain would give rise to a typical Mie plot exhibit- 
ing interference and ripples which are due to the complicated functional form of ajj 
and and not due to the underlying dielectric constant. Surplus electrons would 
moreover not alter the extinction behaviour in this region because |e| 3> |t| and 
|e|»|a|. 

However, immediately above the TO phonon resonance in e, e' < and e" ^ L 
This allows for anomalous optical resonances which are sensitive to small varia- 
tions of e and also to t and a. They are due to resonant excitation of transverse 
surface modes of the sphere II6OI and have been first identified for metallic particles 
where they lie in the ultraviolet |'6r '62l. For a dielectric the TO phonon induces 
them. Fig. 15 shows the dielectric constant e, the refractive index as well as the 
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Fig. 16 Middle panel: Extinction efficiency Qt as a function of the inverse wavelength A^' and the 
radius a for a LiF with = 5 x lO'^cm^^ (left) and an AI2O3 particle with rii, = 3nj./a (right) for 
T = 300K. The dotted lines indicate the extinction maximum for «s = (black), lO'^ (green), 2 x 
10'^ (red), and 5 x lO'^ cm^^ (blue). Top and bottom panel: Extinction efficiency Q, for different 
surface electron densities as a function of A^' and two radii a = 0.2/xm (top) and a = 0.05/Xm 
(bottom). 



extinction Qt, scattering Qs, and absorption efficiency Qa for LiF, CU2O and AI2O3 
particles as a function of A^' and the particle radius a. For LiF and AI2O3 we find a 
clearly resolved series of optical resonances. For the smaller particles they are due 
to absorption while for larger particles they are scattering resonances. The crossover 
between scattering and absorption occurs in the first resonance. For comparison, we 
display in Fig.[T5]also results for CuiO particles which do not show clearly resolved 
resonances. For this material e' < but e" is not sufficiently small to allow for sharp 
resonances. Nevertheless for submicron-sized particles a small extinction resonance 
due to absorption can be identified near 630 cm^' . 

The extinctions resonances are blue-shifted with increasing particle charge f58l. 



This effect is most significant for small particles. Fig. 16 shows the small particle 
tail of the lowest extinction resonance for LiF and AI2O3. The main panel shows 
the extinction efficiency as a function of the particle radius a and the inverse wave- 
length A^^ for «j = 5 X lO'^cm^^ (LiF) or ri/, — 3«i/a (AI2O3). Superimposed is 
the extinction maximum for several values of and corresponding n^. The top and 
bottom panel show the lineshape of the extinction resonance for a — 0.2/im and 




Fig. 17 Position of the extinction resonance depending on the surface charge for LiF, MgO, 
CU2O and AI2O3 (for equivalent bulk charge ni, = 'ins/a) particles with different radii a. Solid 
(dashed) lines are obtained from the Mie contour [Eqs. l |64^ and j65j]. 



a = 0.05/im. For both materials the surplus electrons lead to a blue-shift of the 
resonance. 

For submicron-sized particles where the resonance shift is most significant a < 
X. In this small particle limit we can expand the scattering coefficients for small 
p. To ensure that in the limit of an uncharged surface, that is, for T 0, aJJ and 
b',^ converge to their known small p expansions |54|, we substitute t = T /p prior 
expanding the scattering coefficients. Up to &{p^) this yields a^J = = ^2 ~ ^ 
only b[ ^ i^(p^) contributes. Then the extinction efficiency reads 

e, = I2pie" + a" + 2r'/p) 

{e' + a' + 2- 21" /pf + (e" + a" + 2x'/pf 

where we have restored T. Excess charges enter either by T < 0) or a > 0). 
For T, a — > this gives the limit of Rayleigh scattering. The resonance is located at 
wavenumbers where 



e' + 2-2T'7p=0 for % <0, (64) 

e' + a' + 2 = for j > 0, (65) 

and has a Lorentzian shape provided e" and t' (or a") vary only negligibly near the 
resonance wavelength. Figure [16] confirms the Lorentzian lineshape for AI2O3. For 



LiF, however, e has a hump (see Fig. 15 1 close to the resonance due to a second 



TO phonon mode which is much weaker than the dominant TO phonon. This leads 
to the deviation form the Lorentzian lineshape. 
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For an uncharged surface the resonance is at X^^ for which e' = —2. For X <0 
the shift of the resonance is proportional to t" and thus to n^, provided e' is well 
approximated linearly in X^^ and t" does not vary significantly near Xq^. In this 



case, we substitute in (64i the expansions e' = — 2 + C£(A^' — Aq"') and t" = Citis 
where Cg = de' / {dX^^\^-t and c^- ~ t" /n^l^-t . Then the resonance is located at 

X^^ = Xq^ + Cxn.,/ {TlCgaXQ^). For J > the resonance is located at A^' = Aq"' — 
Canb/ce where Ca = a' /nh\)^-i. 

The proportionality of the resonance shift to Us for LiF and nt for AI2O3 can 



also be seen in Fig. 17 where we plot on the abscissa the shift of the extinction 
resonance arising from the surface electron density (or corresponding bulk electron 
density) given on the ordinate for LiF, MgO {% < 0) as well as CU2O and AI2O3 
(X > 0)- Both bulk and surface electrons lead qualitatively to the same resonance 
shift. Note that even for CU2O which does not show clearly resolved extinction 
resonance a shift is discernible. The most promising candidate for an optical charge 
measurement is AI2O3 where the shift is strongest. 

To illustrate the similarity of bulk and surface electron effects we consider the 
resonance condition or (|65]l for free electrons, which then becomes 



2Ne^ 



3.,2^ =-2 (66) 



e' - 7 ^— ^ = -2 (67) 



for x<Q and 

for X > 0, where is the number of electrons on the sphere. The effect of surface 
electrons is weaker by a factor of 2m* /3m. The factor 2/3 can be understood as 
a geometric factor as only the parallel component of the electric field acts on the 
spherically confined electron gas. Most important, however, the fact that r/p and a 
enter the resonance condition on the same footing shows that the resonance blueshift 
is in the first place an electron density effect on the polarizability of the grain. We 
therefore expect the shift to prevail also for a more complex electron distribution 
on the grain between the two limiting cases of a surface and a homogeneous bulk 
charge. 

Our results suggest to use the shift of the extinction resonance to determine the 
particles charge. This requires the use of particles with a strong TO phonon res- 
onance in the dielectric constant e = e' + /e" which leads to e' < and e" ^ 1 
above the TO phonon frequency. The resoance shift is found for particles with sur- 
face ix < 0, e.g. MgO, LiF) as well as bulk excess electrons (j > e.g. AI2O3). 
For dusty plasmas an optical charge diagnostic can be rather attractive because es- 
tablished methods for measuring the particle charge ll44l |45] '46^ require plasma 
parameters which are not precisely known whereas the charge measurement by Mie 



scattering does not (see Fig. 18 1. Particles showing the resonance shift could be em- 
ployed as minimally invasive electric probes, which collect electrons depending on 
the local plasma environment. Determining their charge from Mie scattering and the 
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Fig. 18 Conventional charge measurement (left) and proposed optical charge measurement by Mie 
scattering (right). Conventional charge measurements rely on trapping the particle in the sheath and 
exploiting a force balance between gravity and the electric force on the particle. The knowledge of 
the plasma parameters is required to infer the particle charge. The optical measurement does not 
require the knowledge of plasma parameters. 



forces acting on them by conventional means B4l l45l l46l would then allow to ex- 
tract the local plasma parameters. Moreover, the Mie signal would provide a charge 
diagnostic for nanodust Ii47l . which is too small and light for traditional charge mea- 
surements. 



5 Summary 

This chapter introduces a microscopic modelling of the plasma wall which com- 
plements traditional sheath models by an interface region — the ESL. In this model 
the negative charge on the plasma boundary is treated as a wall-thermalised elec- 
tron distribution minimising the grand canonical potential and satisfying Poisson's 
equation. Its boundary with the plasma sheath is determined by a force balance be- 
tween the attractive image potential and the repulsive sheath potential and lies in 
front of the crystallographic interface. Depending on the electron affinity that is 
the offset of the conduction band minimum to the potential in front of the surface, 
two scenarios for the wall-bound electrons are realised. For J < (e.g. MgO, LiF) 
electrons do not penetrate into the solid but are trapped in the image states in front 
of the surface where they form a quasi two-dimensional electron gas. For X > ^ 
(e.g. Si02, AI2O3) electrons penetrate into the conduction band where they form an 
extended space charge. 

These different scenarios are also reflected in the electron physisorption at the 
wall. For X < electrons from the plasma cannot penetrate into the solid. They are 
trapped in the image states in front of the surface. The transitions between unbound 
and bound states are due to surface vibrations. In this case the sticking coefficient 
for electrons is relatively small, typically on the order of 10^^. For X > ^ electron 
physisorption takes place in the conduction band. For this case sticking coefficients 
and desorption times have not been calculated yet but in view of the more efficient 
scattering with bulk phonons, responsible for electron energy relaxation in this case, 
we expect them to be larger than for the case of < 0. This indicates that the 
electron affinity should be an important parameter for the charging-up of plasma 
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walls. In particular dusty plasmas offer the possibility to study experimentally the 
charge-up of dust grains depending on this parameter 

Our microscopic model can also be applied to study the optical properties of the 
wall-bound electrons on a dust particle. Surplus electrons affect the polarisability of 
the dust particle by their surface (x < 0) or bulk conductivity (x > 0). This leads 
to a blue-shift of an extinction resonance in the infrared for negatively charged dust 
particles. This effect offers an optical way to measure the particle charge which, 
unlike tradition charge measurements, does not require the plasma parameters at the 
position of the particle. 
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